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Abstract 

Using an elementary example based on two simple harmonic oscillators, we show how a relational 
time may be defined that leads to an approximate Schrodinger dynamics for subsystems, with 
corrections leading to an intrinsic decoherence in the energy eigenstates of the subsystem. 
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It is now well accepted that an energysuperselection rule must be applied to the quantum 
description of the universe as a whole 1]. This implies that the universe does not evolve 
in coordinate time which is usually taken to mean that the quantum state is an energy 
eigenstate. Quite how observers, internal to the universe, develop a notion of time and 
dynamics is then something of a puzzle j2j. One solution is to try and make coherent the 
concept of an internal time based on special subsystems that fulfill the role of clocks. In a 
characteristically elegant presentation Peres js] has given a simple model for how this can be 
done. The Peres model, with three degrees of freedom, reduces to a Hamiltonian which is 
the difference of the Hamiltonians for two simple harmonic oscillators with the constraint 
that the total Hamiltonian vanishes. One oscillator can then be taken as a clock and the 
phase variable for this oscillator is taken as an internal time variable. 

In the example of this paper we consider a Hamiltonian which is the sum of two harmonic 
oscillator Hamiltonians. The constraint is imposed at the quantum level by required all 
physical states to be eigenstates of this Hamiltonian. We show how a relational view of time 
can be given using the fact that such eigenstates are simultaneous near eigenstates of the 
joint phase difference operator. This enables us to treat one oscillator as a clock. However 
quantum fluctuations in the reduced state of the clock oscillator system lead to noise and 
complementary decoherence in the other oscillator. This is made explicit by introducing 
coordinate time using group averaging. A modified Schrodinger equation is then obtained 
for the state of the system which explicitly includes decoherence in the energy basis. The 
approach taken here is similar to the fully relational approach described in |4|. 

While the standard approach takes the state of the universe to be an energy eigenstate, on 
information theoretical grounds, it is more reasonable to assign to the universe a mixed state 
which is diagonal in the energy basis as there is no way any observer inside the universe could 
know which particular eigenvalue is realised. Taking a Bayesian perspective, we need assume 
no more than that the quantum state of the universe is a mixture of energy eigenstates, 



u = 5>|e)<e|. (1) 



P 

At this level of description there is no way to assign the weights p e , so we leave them 
arbitrary. 

An alternative way to write the state in Eq.((TJ) is as a group average^ over the one 
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parameter unitary group generated by the Hamiltonian, 



Pu = lim 1 /V^*)^^ (2) 

where \^f) = ^ e c e |e) is an arbitrary state. In this case p e = \c e \ 2 . In conventional quantum 
mechanics this is the group of time translations and the parameter is coordinate time. 

The universe is a rather large system and admits many decompositions into subsystems, 
some of which correspond to observers. One approach to recovering a notion of time is to 
partition the universe into at least two systems, one of which is to be taken as a clock and the 
rest, which we will simply refer to as the system. It is then possible to decompose the state 
of the universe into a state in which a particular clock variable is seen to be highly correlated 
with particular states of the system. Under the right circumstances this correlation can be 
used to order the states of the system in such a way as to approximate the ordering imposed 
by the one parameter unitary group generated by the Hamiltonian of the system alone. This 
is the so called relational view of time. As we shall see, the time evolution that results from 
a relational approach is not quite the same as unitary Hamiltonian evolution. Indeed, it is 
a special kind of non unitary process, defining a one parameter semigroup, and producing 
an unavoidable decoherence intrinsic to the quantum nature of the universe 0, El . 

As a specific example consider a completely closed system, a 'universe', composed of two 
simple harmonic oscillators. The total Hamiltonian is 



H = uoia\ai + u)2(x\a2 (3) 

(we work in units where h = 1). We have assumed that the oscillators do not interact. 
Anticipating the emergence of a relational time, we now take oscillator labelled 1 as the 
clock and oscillator labelled 2 as the system. 

We now impose a constraint on the allowed physical states of this system: the system is 
in a stationary state of the form Eq.Q where e is the total energy of both oscillators, that 
is to say, we restrict the physical Hilbert space to only those states that are eigenstates of 
the total Hamiltonian. In conventional quantum theory the joint state of two subsystems is 
described by the Hilbert space given by the tensor product of the Hilbert spaces for each 
subsystem. We can use the eigenstates |n)j of the number operators a\ai as a basis for each 
subsystem Hilbert space, 

a\ai\n)i = n\n)i, 71 = 0,1,2,... (4) 



Thus a particular total energy eigenstate can be written as 

|e) = |n)i ® \m) 2 (5) 
with e = c^in + U2m. We now choose uj\ a unit of energy and write the total energy as 

e = n + (D 2 m (6) 

where e = e/coi, cJ 2 = uj 2 /uji. 

For now let us simplify matters by assuming that uj\^ and e are commensurate. In 
particular we assume that 

u 2 = N. (7) 

Together with Eq. (jHJ), this implies that e = M where N, M are integers and n, m are related 
by 

n = M — Nm. (8) 

Clearly the total energy eigenstate is g(M)-fold degenerate with g(M) = + 1- We can 
then write a total energy eigenstate |e) = \M : iV) as 

g(M)-l 

\M : N) = ^ c m \M - Nm)i®\m) 2 (9) 

m=0 

where, as we discuss below, the coefficients c m should be chosen roughly equal to approxi- 
mately recover the conventional time evolution of the system. We note in passing that the 
coherence between the states \M — Nm) 1 ® 1 771)2 superposed in Eq.® is immune to overall 
unitary transformations generated by the total Hamiltonian and is an example of a decoher- 
ence 
Ref. 

We regard the integer M as a label for a particular degeneracy subspace, while the integer 
N is ratio of the system (oscillator-2) frequency to the frequency of the clock (oscillator- 1). 
In other words iV is the parameter that specifies the kind of two oscillator universe we are 
dealing with. The state of this rather sparse universe may then be assigned as the mixture 

Pn = YsPm\M:N)(M:N\ (10) 

M 

where pm is arbitrary at this stage. The assignment of a mixture rather than an energy 
eigenstate is consistent with the idea of a relational universe and follows from the idea that 



ree state for this kind of collective noise, central to the fully relational construction of 
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the particular energy realised by the universe cannot be known. A direct application of 
Bayesian reasoning would then suggests that the most appropriate state of the universe is 
the mixture, Eq. (JTUJ). This of course is quite different to the approach taken in quantising 
ADM canonical gravity in which the Wheeler-de Witt equation implies that the universe is in 
a zero energy eigenstate|2j. However we expect that the more general statistical assumption 
can consistently be incorporated into that approach. 

As above, the one parameter group average Eq. © can be used to relate the coefficients 
c m and pm to an arbitrary state of the two oscillators = J2 nm a n>m |n)i <S> \m) 2 . Indeed, 
we get 



Pm 



g(M)-l 

^ ] a M-Nm,m 
m=0 



CtM-Nm.m /1 1 \ 

and c m = — . (11) 



We now turn to physical quantities that we are permitted to use to describe this universe, 
that is those represented by operators that commute with the total Hamiltonian for the 
double oscillator system as only such operators will be independent of coordinate time. The 
only local operators (that is acting in the individual subsystem Hilbert spaces) that satisfy 
this constraint are functions of the local number operators, a|a«. However when we turn 
to joint operators the situation is more interesting. As we now show the phase difference 
operator is an important example. 

Define the canonical unitary displacement operator for the number basis as 



We write this operator as 



P\n) = \n + 1). (12) 



/•2tt 

P= dfa-^Ety) (13) 
Jo 



where the states E(4>)d(j) is a projection operator valued measure (POVM) defined by 

£(0) = |0>(0l (14) 

with 



|0 = 5>-**|n> (15) 

n=0 

This of course is not a physical state. However we will only ever need the POVM E(<ft) 
which is the optimalgeneralised measurement (POVM) for estimating the phase of a simple 



harmonic oscillator 



il gen< 

tffl. 



Define the operator 

<S= / d(f>E(<j>). (16) 
Jo 

It is then easy to see that 

[$,a + a]=i, (17) 

so we may regard $ and a^"a as canonically conjugate pairs. It is now easy to see that the 
joint operator corresponding to phase difference, $1 — $2 commutes with the number sum 
operator a\ai + a 2 ci2 : a fact which enable one to implement a teleportation protocol based 
on eigenstates of total number 10]. Thus, the operator 

$_ = (T^x - T 2 <f> 2 ) (18) 

with Ti = 2ir/uJi, commutes with the total Hamiltonian in Eq.Q and is accordingly an 
admissible physical quantity given our (self-imposed) constraint. Using the assumed values 
of Di we can write 

$_ = T x ( $1 - -U2 ] (19) 



N 

Now consider the joint phase distribution for the state \M : N), with all the c m chosen 
as equal, c m = g^ 1 ^ 2 

P(<j>xM = tr (Exfa) ® E 2 {ct> 2 )\M : N)(M : iV|) (20) 
Using equations Eq. p4jl and Eq.(|15jl we find that 



P(0i,0 2 ) = ~ 
9 
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1 _ e iWi-W 

When M > iV, or equivalently ^ 3> 1, this is sharply peaked at 



(21) 



JV0i - 02 = 2vrA;, k integer. (22) 

Given the definition of $_ in Eq.()19|l. we see that the state \M : N) is a near eigenstate of 
the phase difference operator. 

The interpretation of Eq. (J22j) goes as follows. Independent measurements of the phase 
of each oscillator will reveal a strong correlation. Given an arbitrary "clock reading" <fii, 
the conditional probability of the system's phase 2 is sharply peaked around the value 
N(pi mod 27T. This is a very familiar structure that can be more easily seen using a graphical 
construction shown in figure [T] a). We end up with the Poincare section in the phase space 



of the two oscillators for which 02 is a function of 0i. Classically this immediately implies a 
dynamical evolution with respect to a coordinate time t = <pi/ui that is precisely what one 
would expect. 




FIG. 1: a) A schematic representation of the phase space solution for 0i and 02 in the asymptotic 
limit g — » oo, with N = 3. b) The joint distribution of 0i and 02 given by Eq. ()20|) for the 
parameters N = 3 and g = 14. As expected, the distribution is sharply peaked around 02 = 
N(j)\ mod 2tt, so the system undergoes three complete oscillation for 0i £ [0, 2tt]. These correlations 
inherent in the phases of the two oscillators in a near eigenstate of phase difference. As M increases, 
the distribution reaches a delta function and the relation 02 = ujit mod 2tt becomes exact [as 
illustrated in figure a)], hence the definition t = (pi/uoi. 

However as M is in fact finite, the correlations between 0i and 02 are imperfect, as seen 
on figure Qb). When viewed as a dynamical system and insisting on recovering differential 
equations that are local with respect to the "time" variable 0i, these fluctuations will be 
interpreted as noise or error in the frequency of the system. In other words, quantum 
fluctuations inherent in the nature of the total state of both oscillators \N : M) will appear 
as an intrinsic source of noise in the classical record of measured phases. This implies an 
effective decoherence in the local energy eigenstates for the system, as we now explain. 

Returning to the equivalent expression for pu given in Eq.(j2J) , but now written as a 
discrete average 

1 n 

pu = d 1 ™ rTTi E e " m ^>^l ® l£>2<£l^ (23) 

we can give a presentation of relational time in terms of a coordinate time t. In this picture it 
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is natural to interpret each term in Eq. (j23J) as the coordinate time evolution with 1^)1 ® |£)2 
as an initial state of the two oscillators when they do not interact. We stress that these 
initial states are completely arbitrary. However from a perspective internal to the universe 
in which we wish to interpret oscilllator-1 as a clock, we may only wish to consider those 
situations in which \ip)i is chosen especially. How is this to be done? 

A good clock must have some dynamics, in fact we would like it to evolve rapidly compared 
to everything else. This means it must be in a state with a large variance in energy, which 
here implies a large variance in number a\ai. So we will take \ip) so that 

(^{alatfm - ((Ma^y » 1 (24) 

which means that c n = i(n\ip)i is almost independent of n over some large range K <n < 
K + L with L ^> 1. In this case, everything else means simply oscillator-2. We do not specify 
the state of this system at all. Our objective is to recover standard Schrodinger dynamics 
for this system. 

Suppose we now ask for the conditional state of oscillator-2 selected according to all those 
measurements of $ on the clock system that gave a particular result <p\. This is given by 

= tn{\ct>)M®i2Pu) (25) 

Substituting the expression in Eq. (j23j) and using H = a\a\ + Na\a 2 we find that 

n 

pf> = V{<f)\a)e- iaNa ^ \02(£\e taNa * a2 (26) 



where 



(27) 



2 

V{(j)\a) = W)]' 1 ^ ' »-« a -^ n 

n=0 

where the normalisation is the inverse of the probability that the result is indeed selected 

n 

^) = (HlEW (28) 

a=0 

The form given in Eq.(|26|) is the same construction as recently proposed by Gambini et 
al. in their construction of relational time[6]. It implies the intrinsic decoherence that we 
alluded to the the relational construction of the Poincare section. 

When is chosen to be a good clock state according to the prescription Eq. (j24|) we 
can show that V(4>i\a) considered as a function of the integer a is sharply peaked at a = 4>i- 

8 



In this case we may approximate the selected state as 

pf^ « e -*Wi<4<w|£) 2 (£| e *Wi<4<w ( 29 ) 

using a simple redefinition of the variable <pi = U\t in terms of a 'time' coordinate and using 
the fact that 002/^1 — N we see that this equation may be written 

pf l] = p 2 {t) w e - iLJ2ta 2 a2 |02(^|e i " 2ia 2 a2 (30) 

which of course is Schrodinger dynamics with respect to coordinate time t. Successive 
corrections to this approximation describe an intrinsic decoherence in the energy basis. 

To make further progress we need to specify a clock state and compute V(<f)\a). How- 
ever a particularly simple model that replaces unitary Schrodinger evolution with a unital 
semigroup evolution was first proposed in [7] . This is done by taking a Poisson distribution, 

pU\a) = M! e -^ (31) 
a\ 

where 7 is dimensionless number which parameterises different clock states. Making the 
change of variables, uj\t = (ft, 7 = 7/Vl we can write this as 

V(t\a) = M^e" 7 * (32) 

where now 7 may be taken as a new fundamental constant with dimensions of frequ ency 
that takes a maximum value and limits us in principle to a best, but imperfect, clock 12]. 
In the presentation given above, 7 determines how fast the dynamics of a clock can be by 
restricting the allowed variance in the energy of the clock to some finite maximum value. 
One expects that this limit will arise naturally in a fully quantum theory of the universe, 
that is, one in which both matter and spacetime are given a quantum description. However 
at this stage it is left as a free parameter. Of course we might expect (7) -1 to be related 
somehow to the Planck time. 

Using Eq.(jH2J) we can easily see that the Schrodinger equation is replaced by 

= 7 (e-^/>(t)e*4« 2 /7 - p 2 (t)) (33) 

where we see that the right hand side is the generator of a unital semigroup that leaves 
all quantities that commute with the Hamiltonian invariant in coordinate time. If, as we 
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expect, 7 is very large, but not infinite , we can expand in inverse powers of 7 to obtain a 
first order correction to Schorddinger dynamics 

= -*[4 a 2,P2(*)] - 7^[4 G 2, [4«2, P2(t)}} + ■■■ (34) 



The second term here leads to diffusion in the phase variable $2, as suggested in the 
relational Poincare section, while coherence in the number basis is damped. The general im- 
plications of this dynamics are explored in Q. One important consequence is a modification 
of the dispersion relations for field theories. 

In our construction the emergence of coordinate time, and an intrinsic decoherence pro- 
cess, arises due to the freedom to choose states of the total system as mixtures of energy 
eigenstates. This freedom incorporates the same intuition that suggests that the state of 
the universe must be invariant under the quantum time translation operator. While the 
universe must be in an energy eigenstate, it is seems likely that the particular eigenvalue 
realised cannot be known. A Bayesian perspective would then suggest that we assign the 
state of the universe as a mixture over energy eigenstates, with an arbitrary distribution. It 
is this property that enables us to introduce coordinate time as a group average. It would be 
interesting to consider groups other than the time translation group, for example the Lorentz 
group. The fact that the mixture is arbitrary gives us the freedom to describe subsystem 
dynamics in coordinate time from arbitrary initial conditions. The choice of a mixture of 
energy eigenstates is a departure from what is usually done in canonical quantum gravity. It 
would be interesting to revisit the assumptions in canonical quantum gravity to understand 
this freedom. 
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